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ZrZn2: geometrical enhancement of the local DOS and quantum design of magnetic 
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The recent discovery of coexisting ferromagnetism and superconductivity in ZrZn^, and the fact 
that they are simuhaneously suppressed on applying pressure |l| suggest the possibihty of a pairing 
mechanism which is mediated by exchange interactions and connected with the proximity to a 
magnetic quantum critical point. On the basis of first principles, full potential electronic structure 
calculations, we study the conditions that, for ZrZn2, determine the proximity to this magnetic 
instability. More specifically, we discuss the role played by the geometrical arrangement of the 
lattice, hybridization effects and the presence of disorder, as well as the application of external 
pressure. These circumstances influence the width of the relevant Zr d bands whose narrowing, due 
to the reduction of the effective number of neighbours or to an increase of the cell volume, causes an 
enhancement of the density of states at the Fermi level. Finally, we highlight some general features 
that may aid the design of other materials close to magnetic instabilities. 
PACS:74.25.-q,74.25.Ha,71.20.Be 



The very recent discovery of the coexistence of ferro- 
niagnetism (FM) and superconductivity (SC) in ZrZn2 
and UGe2 suggests that the current ideas on SC 
need to be critically reconsidered. Interestingly ZrZn^ is 
a relatively 'simple' metal, in the sense that there are not 
complications arising from / electrons |3|] or oxide planes. 
This circumstance implies that its electronic properties 
can be studied in detail within the density functional 
theory M framework which is expected to be reliable for 
such systems. 

As known for many years, ZrZn^ is a weak itinerant 
FM [g| with a Curie temperature TpM = 28.5if , despite 
the fact that both pure Zr and Zn are not magnetic. 
Both pure elements also exhibit superconducting transi- 
tions with Tc = 0.6K and Tc = 0.85K, respectively [||. 
Many speculations on the properties of this compound 
have been stimulated by the discovery R that, at some 
critical pressure, Pc — SKBar, it becomes paramagnetic 
(PM), suggesting that it is very close to a magnetic quan- 
tum critical point (MQCP). Notwithstanding the argu- 
ment of Berk and Schrieffer ||8|] about the suppression 
of the phonon-induced s-wave SC by FM spin fluctua- 
tions, many authors ||9|] argued that, in the vicinity of 
the MQCP on the PM side (high pressures), p-wave SC 
could be induced by paramagnons. Later on, Fay and 
Appel ||I^ suggested the possibility of having p-wave SC 
also for P < Pc , inside the FM region of the phase dia- 
gram. A different scenario, recently proposed by Blagoev 
et al. uM, predicts s-wave SC for P < Pc and p-wave SC 
for P > Pc- 

The recent experiments M suggest that the kind of 
SC observed in ZrZn2 is related to the presence of FM, 
since both FM and SC are suppressed simultaneously as 
the pressure reaches Pc- This new phenomenon could 
be compatible either with the picture of Rcf. flo] or with 
that of Rcf. [0 , both assuming a magnetic excitation- 



mediated coupling in the proximity of the MQCP. In 
order to decide between alternative explanations it is 
highly desirable to be able to engineer materials close 
to a MQCP and to control the vicinity by changing some 
variable such as the external pressure or the concentra- 
tion of some impurities. In this letter, on the basis of 
first-principles, full-potential, electronic structure calcu- 
lations for both ZrZn2 and related systems, we highlight 
a rather general feature which may aid this endeavor. 

In particular we examine the ZrZn2 lattice geometry 
and follow some important consequences for the physical 
properties. We find that the density of states (DOS) 
at Fermi energy, n{Ep) is enhanced, in a sense to be 
specified below, by the lattice geometry. A high n^Ep) 
is a crucial attribute for a material to be close to an 
MQCP as shown, for example, by a simple Stoner factor 
calculation (a measure of proximity to a ferromagnetic 
phase transition) 



So^{l^In{EF)y^ 



(1) 



where / is the exchange integral. Of course, a large n{Ep) 
is also important in generating a sizeable electron-phonon 
coupling, 



A 



njEp) <P> 

M <LU^ > 



(2) 



where < /^ > is the electronic stiffness parameter of 
Gaspari and Gyorffy [Q, M the ionic mass and < uj^ > 
some mean value of the phonon frequency. In the above 
equations we have deliberately used definitions from the 
simplest available electronic theories for magnetism and 
SC, just to stress the relevant role played in these phe- 
nomena by n(Ep). 

ZrZn2 crystallizes into a cubic C15 superlattice p3J. 
In this structure, Zr atoms occupy the positions of a di- 



amond lattice while the Zn atoms form a network of in- 
terconnected tetrahedra. Since the major contributions 
to n{Ep), as we shall see, come from Zr, the local envi- 
ronment of Zr atoms (Fig. [l|) is particularly important 
for our concerns. Each Zr is surrounded by 12 Zn neigh- 
bours and 4 Zr neighbours, at distances that in terms of 
the lattice constant, a, are, respectively, dzrZn — "g"'^' 

and dzrZr = 40- In other words, the distance between 
two Zr atoms is only 4 per cent larger than the near- 
est neighbours (NN) distance dzrZn- This circumstance 
leads to an appreciable overlap of the wave-functions of 
two neighbouring Zr atoms |14[| . 

We have calculated the band structure of ZrZn2 in the 
PM state using the full potential linearised augmented 
plane waves method [Ea (FLAPW) within the local den- 
sity approximation (LDA). The calculated equilibrium 
lattice constant (13.58 a.u.) is about 3 percent smaller 
than the experimental one (13.98 a.u. ||]). Importantly, 
our full potential results provide a reference for more ap- 
proximate, existing band structure calculations |l4,E6[, 
all based on the linearized muffin-tin orbital (LMTO) 
method and the atomic sphere approximation (ASA). 
These simpler band structure calculations form the start- 
ing point for studies of electron-phonon coupling and spin 
fluctuations in this material. 

The DOS (Fig. 0) classifies this compound as a split 
band system: the electronic states are distributed into 
largely separate energy regions. The Zr d resonances are 
close to Ep, while those of Zn occur well below. The band 
structure, plotted in Fig. pi reveals very narrow d bands 
and suggests that these would be appreciably broadened 
out by a reduction of the interatomic distances. A sim- 
ilar broadening could also be caused by decrease of the 
quasiparticle lifetimes due to impurities or defects. This 
view is supported by Fig. H, where we see how the in- 
crease of the bandwidth for smaller values of a results 
in a dramatic drop of the DOS peaks heights. For the 
system at hand, the drop of n{EF) is also affected by elec- 
tronic topological transitions |17| (ETT) occurring when 
the Fermi level crosses the various DOS peaks, as shown 
in Fig. ^. 

The variation vs. a of the local DOS in Zr muffin- 
tin spheres at the Fermi level, nzriEp), is displayed in 
Fig. ||. These trends can be easily fitted using a simple 
tight-binding (TB) model [Q: nzr{Ep) is inversely pro- 
portional to the bandwidth which is proportional to the 
hopping parameters. These scale exponentially with the 
relevant distances, say as e~^^ with _B > 0. Accordingly 
with Eq. |l|, the behaviour of nzriEp) is consistent with 
the experimental findings 0M] that a compression can 
suppress the magnetic ordering by crossing, at P = Pc, 
the MQCP (5(7^ = 0). On the other hand, consistently 
with Ref. ||l|, the MQCP could be reached at lower pres- 
sures in impure samples, due to the lifetimes effects men- 
tioned above. 

Pure Zr forms a hep crystal structure. Although 
very well known, the standard explanation as to why 



it is not ferromagnetic contains some valuable pointers. 
nzr{Ep) = 4.11 is about 4 times smaller than in ZrZn2. 
Such a large factor cannot be simply attributed to the 
atomic distances involved (see Fig. g). For pure (hep) 
Zr wc find dzrZr — 5.900 a.u., very close to both dzrZr 
and dzrZn in ZrZn2- Neither can this factor be caused 
entirely by ETT's. Such a decrease can be explained only 
by assuming that, as a consequence of the very weak hy- 
bridization with Zn, Zr atoms 'see' only their 4 Zr neigh- 
bours, while in pure Zr the number of NN's, /, is 12. 
This could account for the reduction of the local DOS by 
factor 3. Using the above TB model, we can draw up a 
relationship 



nzriEp) - Ae'"'-'-''/f,ff, 



(3) 



where f^ff is some effective number of NN's. Of course, 
Eq. y should not be taken too literally since ETT's and 
differences between the shapes of the DOS's in different 
lattices (e.g. the pseudogap between bonding and anti- 
bonding states in hep lattices) have also some influence. 

In order to add substance to our picture of ZrZn^ in 
which the Zn atoms play essentially the role of empty 
spheres placed between Zr atoms to reduce j^f j , we have 
calculated the bandstructure of Zr in the hypothetical di- 
amond lattice phase that would result by the substitution 
of the Zn atoms in the ZrZn2 lattice with empty spaces. 
The comparison of the corresponding DOS (Fig. m) shows 
up some minor differences mostly due to the weak hy- 
bridisation with Zn and with the charge transfer to Zn 
atoms in the ZrZn2 system. 

In summary: we have pointed out how the addition 
of Zn to an ideal Zr diamond lattice changes the local 
DOS on Zr atoms very little. We surmise that this is an 
example of a rather common occurrence. One can con- 
sider a hypothetical geometrical lattice structure formed 
by adding empty spaces in some metal. The narrowed 
relevant bands could give rise to exotic, magnetic or su- 
perconducting, properties (as shown here for Zr in the 
diamond lattice, that, of course, cannot be realized in 
practice) . Those same properties could then be obtained 
by filling the empty spaces with atoms which give rise 
to little (or zero) DOS at the Fermi level (as Zn) and 
which have little hybridization with the metal we want 
to modify. Contrary to the fashion prevalent in much sci- 
entific literature, we close this Letter by saying that the 
observations we are making from our calculations are not 
particularly new. The focus, however, on the geometrical 
guideline we have extracted is. A list of existing systems 
whose properties seem to be consistent with this guideline 
can be easily constructed. For example, there are thin 4d 
and 5d metal films ||l^, where the reduction of /e// at 
the surface allows magnetically ordered states to be sta- 
ble. Another example might be the alkali fuUerides [ pO[ , 
where Cgo molecules could provide the 'vacuum' mak- 
ing alkali metals superconducting. Other possible exam- 
ples are MgSa Hi, f/Gea §, MnSi 0,|| the layered 
cuprate p^ and ruthenate superconductors ||23] and A\h 
compounds M. 
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FIG. 1. The local environment of a Zr atom (marked in 
black) in the C15 cubic Laves structure. Other Zr atoms: 
dark circles; Zn: red circles 
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FIG. 2. DOS of PM ZrZn2 (top panel) at the relaxed lat- 
tice constant, a=13.58 a.u.. The labels Zn and Zr refer to the 
corresponding mufRn-tin spheres. The lower panel shows the 
interstitial DOS divided by the number of atoms in the unit 
ceU. 
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FIG. 5. DOS of PM ZrZni about Ef for the values of a 
indicated by the labels. 



W L r XWK rXUWKLU 

FIG. 3. Band structure of PM ZrZn^ along relevant sym- 
metry directions. 





FIG. 6. n{Ep) vs. dzrZr in ZrZn2. The circles correspond 
to the calculated values, the line is an exponential fit. We 
find B=1.2 (see the text), in good agreement with the values 
reported in Ref.[18]. 
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FIG. 4. DOS of PM ZrZn2 for the values of a indicated 
by the labels 
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FIG. 7. The density of states of pure Zr in the diamond 
lattice (dots) compared with the density of states of ZrZn-i 
(plain line), at the relaxed lattice constant of ZrZn^. All the 
DOS refer to Zr muffin tin spheres only. The arrow indicates 
the direction in which the Fermi level is espected to move due 
to charge transfer from Zr to Zn. The dashed line marks an 
upper bound estimate of the above Fermi level shift. 



